Read out of a Nuclear Spin Qubit 
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We propose a detector to read out the state of a single nuclear spin, with potential application in 
future scalable NMR quantum computers. It is based on a "spin valve" between bulk nuclear spin 
systems that is highly sensitive to the state of the measured spin. We suggest a concrete realization 
of that detector in a Silicon lattice. Transport of spin through the proposed spin valve is analogous 
to that of charge through an electronic nanostructure, but exhibits distinctive new features. 
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Nuclear magnetic resonance (NMR) experiments have 
been a valuable testbed for quantum information process- 
ing (QIP) and they still provide the largest collections 
of coupled qubits available at present 0, • Most NMR 
QIP experiments are performed on liquids and suffer from 
the lack of scalability. Solid state spin systems have been 
roposed as a promising route to scalability in NMR QIP 
Their experimental implementation is, however, 
challenging. One major obstacle that has to be overcome 
is the read out problem. In most proposals it requires the 
measurement of the quantum state of a single nuclear 
spin. Experiments that have successfully detected single 
spin resonances are promising [3, H Oj floL . Recently 
the read out of a single electron spin has been reported 
[l2j| . The measurement of the state of a single nuclear 
spin has remained, however, elusive up to now. The adi- 
abatic transfer of the spin state of a nucleus to that of an 
electron or magnetic resonance force microscopy 0] 
have been proposed to read out an NMR qubit. These 
techniques introducejhowever, unwanted sources of addi- 
tional decoherence [14|. The use of ensembles of nuclear 
spins in 29 Si as qubits has been proposed to enhance 
the measurement signal. Optical detection by measuring 
the energy of photons emitted by bound excitons is an- 
other proposal to measure a single nuclear spin 0] . Most 
of the above mentioned schemes require very specific ma- 
terial properties for their implementation. In this Letter 
we propose a scheme for the read out of a nuclear spin 
qubit that relies only on the dipolar coupling between 
nuclear spins. That interaction is generically present in 
solids. 

Our proposal is inspired by charge detectors that uti- 
lize quantum point contacts (QPC) [H E [HE 13 
l2(il ]. Such a detector has recently been successfully em- 
ployed for the read out of an electron spin qubit after 
spin to charge conversion |l2j | . Converting a nuclear spin 
into a charge signal is hard and we therefore propose to 
directly measure a nuclear spin by an analogue of a QPC 
for nuclear spin currents, as shown in Fig. ^ Our de- 
tector consists of two bulk systems of nuclear spins sf^ R 
connected by two spins Sj that act as a spin valve. A 
difference in polarization of the bulk spins to the two 
sides of the valve can drive an equilibrating spin current 
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FIG. 1: Proposed nuclear spin detector. Upper panel: 
schematic setup with inter-spin coupling constants. Lower 
panel: Possible realization in a Silicon lattice. Filled circles: 
29 Si. Open circles: 28 Si. 



between them. A qubit spin I q and an auxiliary spin l a 
(that is prepared in a known state) create local magnetic 
fields for Si and S2. These fields depend on the qubit 's 
state and can bring the valve spins into and out of reso- 
nance. We demonstrate below that as a consequence the 
spin current between the bulk systems is highly sensitive 
to the state of the qubit. After an appropriate measure- 
ment time the state of the qubit is therefore encoded in 
the spin state of a measurable number (typically 10 6 ) of 
bulk spins and can be read out with standard techniques. 

We write the Hamiltonian of the spin detector Fig. 
as the sum 



H — Hb + H v + Hq 

L/R 



(1) 



We assume the bulk spins s^- to form a linear chain of 
sites (N is large) with Hamiltonian 
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sf*sf* 1+S f y s% 1 + {L~R} 



(2) 



H v describes the valve spins Sj coupled to the qubit l q , 
H v = 2H{S?S$ + SfS y 2 ) + B(I*S1 + I*S Z 2 ), (3) 

and Ht couples the valve spins to the bulk, 

H T = 2t[s^(Sf+aSI) + s^(Sf + aS^)} 

+ 2t[ S K*(SZ + aS?) + S K y (S% + aSV)]. (4) 
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The choice of coupling of the qubit to the spin valve 
BI*Sf is crucial. It assures that the qubit's state is con- 
served after having been projected onto the logical basis 
{It), I 4)} (^IT) = I T>/2 and /|||) = -|4)/2). This 
allows the detector to be operated until it has accumu- 
lated a detectable change in bulk polarization. Experi- 
mentally this form of the coupling Hamiltonian can be 
implemented by choosing the detector nuclei to be of a 
sort different from that of the qubit nucleus. A difference 
of the gyromagnetic ratios of the two kinds of nuclei to- 
gether with a strong magnetic field along the z-direction 
then eliminates the couplings I*Sf and I]fSf that are 
otherwise present. It removes the corresponding transi- 
tions far off resonance and leads to an effective Hamil- 
tonian of the form H v . By means of a Jordan- Wigner 
transformation the spin operators in H can be expressed 

L / R 

in terms of fermion operators c- and dj with stan- 
dard anti-commutation relations |2l| . In Fourier space, 
c r '(p) = (l/V^EjcJeaptpj (77 € {L,R}), the Hamil- 
tonian takes the form 

H b = b co8pd^(p)tf>(p) 
v e{L,R}, P 

H v = n&d 2 + h.c.) + B[%{£d 1 -\) + IZ&d 2 -\)\ 

H t = J2 + (5) 

ve{L,R} 

Both I q and I a are not dynamical and H b and H v are 
readily diagonalized in their fermionic form. The non- 
linear terms 

Tjj = t(d 1 +ae l ^d 2 ), 

T R = t(d 2 +ae~ M ^d 1 ), (6) 

however, introduce an interaction between the fermions. 
This precludes a straightforward analytic solution of the 
transport problem. As we show below, these interac- 
tion terms affect the transport behavior of the spin valve 
in a qualitative way. In Ref. |22j | similar interactions 
have been encountered for hard core bosons. There they 
were dealt with numerically. Here we exploit the fact, 
that in typical NMR experiments the temperature T is 
high compared to all intrinsic energy scales. This allows 
for a modified mean field treatment of the bulk spins 
that can be carried out analytically. In this limit, the 
bulk spins rearrange to their equilibrium state instantly 
after every transfer of spin into the bulk. To a good 
approximation their dynamics is therefore independent 
of that of the valve spins. This reduces the complex- 
ity of the problem to the four dimensional Hilbert space 
of the two valve spins Sj. It leads to a quantum mas- 
ter equation for the reduced density matrix of the valve 
spins. Such an approach often proves useful in problems 
of quantum transport |24j . It is valid if the time scale 
t c = 1/ max{f2, B, t 2 /b} of the dynamics of Sj is much 
longer than that of the bulk dynamics r b = 1/ min{T, b}, 

T c > T b . 

In deriving this master equation we proceed closely 
along the lines of Ref. We only repeat the main 

steps here. We evaluate the generating function 

Z(\) = Tr e - 4A ^" e - lHT e iA " s " p (in) e lHT . (7) 



p( m ) is the initial density matrix of the spin system and 
gL/R _ s L .l Rz i s the total bulk spin. Z generates 
moments of the amount AS M of spin that is transferred 
through the valve spins during time r. In particular, the 
mean amount of transferred spin averaged over identical 
experiments is 

(ASr>)=i-£- lnZ(A) . (8) 

dXrf A=0 

It is crucial to note, that the Hamiltonians H in Eq. (J7J) 
as well as S v are quadratic in the bulk fermion operators 
c^. Hence they can be integrated out exactly, yielding 

Z(\) = Ty c T ± p { i n) e~ l S * dtdt ' S„ €{£) R } rt(*)Gf„(t-f)f,(t') 

x e ~i Jo dt n{d\T*d 2 +h.c.)+B(I* d\ T *d x +Il 4r'd a ) _ /g\ 

Here, all operators are vectors in a "Keldysh space" 
(dj, dj) of operators dj and dj that originate from the 
first and the second exponential of Eq. (7J respectively. 
By the symbol T± they are ordered in time as well as 
relative to the initial reduced density matrix pc of the 
central spins Sj and I 9 /I Q . r 2 is the third Pauli matrix 
and Tjj = t z exp(— iX v T z )T v . The "mean field" due to 
the bulk spins is able to increase or decrease Sj and it is 
quantified by the bulk spins' Green functions G v . G I? (i) 
are peaked around t = with width 17,. In the high tem- 
perature regime r c 3> t of interest here the exponent in 
Eq. © is therefore local in time and Z is the integral of 
an ordinary differential equation. We write it accordingly 
as the trace over a time-dependent density matrix, 

Z(X) = Tr c P x c (t), (10) 
that obeys the master equation 
d tP x c = C[p x ] = -1 [H v+Ptl (TX + T 7] T^),p x c ] 

~ Y h l n v T v T nPc + C 1 - n v )T^T v p x + h.c] 
+ 1 [n v T}p x T n e lA " + (1 — n v )T v p x T^ e^] 

(11) 

with initial condition p x (0) — p™ . Due to the coupling 
to the bulk spins, p x is evolved in time by a linear "su- 
peroperator" C of Lindbladian form [25} rather than by 
a Hamiltonian. n v — 1/2 = Trp ln SQ Z are the expecta- 
tion values of the z-projections of the bulk spins and 
Prf = (2Trb)^ 1 V J dujn^/uj are principal value integrals 
that are cut-off at high frequencies by the bandwidth b 
of the bulk spin excitations. We compute Z at large 
times t by exponentiating the largest eigenvalue of tC. 
The spin currents j n — AS n /t then follow from Eq. JSJ. 
We take an initial density matrix p™ corresponding to I a 
being prepared in state |f). If the qubit's state is ||), the 
spin current has magnitude 
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jf = 2(1 + a 2 )Anr 



a 2 B 2 + (1 - a 2 ) 2 [a 2 r 2 (l - £An 2 ) - 2£afi£An(l - a 2 



(1 + a 2 ) 2 B 2 + (1 - a 2 ) 2 [(l + a 2 ) 2 r 2 - 4(afi5An(l 
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(12) 
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FIG. 3: Transport processes for a typical initial spin config- 
uration to lowest order in perturbation theory. Large circles: 
spins Si and S^. Small circles: bulk spins So and s$. Filled 
circles symbolize spin up (occupied site), empty circles spin 
down (unoccupied site) for a spin (fermion) system. 



FIG. 2: Measurement contrast C of our detector as a function 
of the strength B of its coupling to the qubit - normalized to 
the coupling Bo to a proton qubit (solid line). Dashed line: 
C for a corresponding system of free fermions. 



Here, Y = t 2 /b, Q = fl + Aat 2 {p L + p R ), An = n L - 
nji, and £ = 1. If on the other hand is in state ||), 
the magnitude of the spin current is jf = jf | B=0 . The 

measurement contrast C, denned as the ratio C = jf/jf 
of the signals for both possible states of I q , characterizes 
the performance of the detector [2a] . 

Fig-fflsuggests an implementation of our proposal in a 
chain of Silicon atoms [27]. The weak links between the 
bulk and the valve spins are realized by lattice vacancies 
or isotopes 2 8Si without nuclear spin. In typical solids, 
nuclear spins are coupled by the dipolar interaction that 
falls off as r~ 3 with the distance r between spins. For sin- 
gle vacancies we therefore have t = Q = 6/8 and a = 1/8. 
In Fig. [21 we plot the measurement contrast C for these 
parameter values at full polarization (rij, = 1, n R = 0) as 
a function of the coupling strength B. B is normalized to 
the value Bq that one has for the measurement of a pro- 
ton spin. We estimate B = (8/5) 3 (6/4). Fig. clearly 
demonstrates that our detector yields good contrast over 
a large range of coupling strengths. In particular, we 
find excellent contrast C w 60 for the measurement of a 
proton spin. 

For full polarization the measurement contrast can in 
fact be made arbitrarily large by fine tuning B. This 
important property of our detector can be traced to the 
presence of the interaction terms in Eqs. ©. This fol- 
lows from a comparison of C for our spin system with 
that of the free electron system described by H with- 
out the phases exp(i7refid») in T v (that descibes a double 
quantum dot |24|)- C in that case is obtained from Eq. 
(|12fl by setting £ = and it is bounded from above, as 



shown in Fig. [21 The discussed divergence of C for the 
spin valve is due to an interference effect that causes 
to vanish at 



B mt - (1-a 2 )-. 

a 



(13) 



At this coupling strength, there occurs a complete de- 
structive interference of transport processes involving 
next-nearest neighbor couplings at with processes involv- 
ing nearest neighbor couplings t and Q only. To illustrate 
this, we show in Fig.[3all lowest order transport processes 
for a typical initial spin configuration. Adding up their 
amplitudes leads to a transport rate cx \t 2 (fl/B — a)\ 2 
that indeed vanishes at B = B int (to lowest order in a 
and t). For free fermions the same processes exist with 
the correspondences spin-up <-*■ occupied site, spin-down 
<-> unoccupied site. However, the last process d) acquires 
an additional minus sign, because the fermions in the fi- 
nal state are interchanged with respect to those in the 
final state of process b). This results in a strictly non- 
vanishing rate cx \t 2 fl/B\ 2 that is symmetric under the 
reversal of the sign of B. These features carry over to 
all orders of perturbation theory and render transport 
through the spin valve qualitatively different from that 
through the corresponding system of free fermions, as 
seen clearly in Fig. [21 

In case the polarization of the bulk spins is not com- 
plete, C does not diverge anymore. Instead, it possesses 
a maximum at an optimal coupling strength B max . We 
plot C m ax and B max as functions of the polarization P 
in Fig. E| taking n L = (1 + P)/2 and n R = (1 - P)/2. 
Evidently good contrast can be attained even for a weak 
polarization of the bulk spins. Note, though, that the 
plotted maximal contrast is in practice increasingly hard 
to achieve for decreasing polarization, since the necessary 
optimal coupling B max diverges at P = 0. 

We finally comment on the applicability of our model 
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FIG. 4: Maximal contrast C max (solid line) together with the 
the optimal coupling strength B max at which it is attained 
(dashed line) as a function of bulk polarization P. 

to real spin systems. We have analyzed a one-dimensional 
spin chain and have neglected next-nearest neighbor cou- 
plings as well as s 1 j Z s^ 1 interaction terms for the bulk 
spin systems. In principle, effectively one-dimensional 
spin chains are available in fluorapatite |28j and proposed 
in Silicon by Itoh. The mentioned couplings can be elimi- 
nated with experimental effort |2^|. We believe, however, 
that our predictions remain qualitatively correct also in 
more general situations. Spin systems with sj z in- 
teraction are described by interacting fermions [3(J]. To 
lowest order in t our Eq. 10 holds also in the presence 



of bulk fermion interactions, with being the Green 
function of the interacting system. For one-dimensional 
interacting fermions G v exhibits an anomalous behavior 
that is cut-off at low energies by the temperature of the 
system. In our limit of temperatures that are larger than 
the bandwidth we do therefore not expect that these in- 
teractions have a qualitative effect. The same reasoning 
applies to possible next-nearest neighbor couplings and 
higher dimensionality, that are similarly described by ad- 
ditional interactions between bulk fermions in our model. 
Inclusion of the direct dipolar coupling between the spins 
Sq and s§ leads to a leakage current that is suppressed by 
a factor of 3~ 6 10~ 3 relative to the leading contribu- 
tion. The coupling S* Sf has no effect on spin transport 
and to lowest order in t the couplings s rlz Sj result in a 
trivial shift of the energies of S?. 

In conclusion, we have proposed a local read out 
scheme for future solid state NMR quantum computers. 
It is based on a spin valve between bulk nuclear spin 
systems. Our analytical results show that it is highly 
sensitive to the state of a nuclear spin qubit. It imple- 
ments a detector of single nuclear spins. An experimental 
realization of the proposed detector would not only solve 
a major problem on the way to a scalable NMR quantum 
computer. Our comparison with the corresponding elec- 
tronic system indicates that it would also provide us with 
a system that exhibits novel and interesting phenomena 
in quantum transport. 
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